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1 Abstract. The fire history of the northern larch forests within the permafrost zone in a portion of northern 

2 Siberia (~66°N, 1 00°E) was studied. Since there is little to no human activities in this area fires within the 

3 study area were mostly caused by lightning. Fire return intervals (FR1) were estimated based on bum 

4 marks on tree stems and dates of tree natality. FR1 values varied from 130 yr to 350 yr with 200±50 yr 

5 mean. In southerly larch dominated communities FRi was found to be shorter (77±20 yr at ~61°N, and 

6 82±7 at 64°N), and longer at the northern boundary (~71°) of larch stands (320 ±50 yr). During the Fittle 

7 lee Age period in the 16th to 18th centuries FRI was approximately twice as long as recorded in this 

8 study. Fire caused changes in the soil including increases in soil drainage and permafrost thawing depth 

9 and a radial growth increase of about 2 times (with more than 6 times observed). This effect may simulate 

10 the predicted wamiing impact on the larch growth in the permafrost zone. 

11 
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1 1. Introduction 

2 Larch ( Larix spp.) forests compose about 43% of Russian forests. Larch stands are dominants from the 

3 Yenisei ridge on the west (~92°E longitude) to the Pacific Ocean and from Baikal Lake to the south to 

4 73rd parallel to the north. Wildfires are typical for this area with the majority occurring as ground fires 

5 due to low forest crown closure. Larch is a pyrophytic species; since fires promote the establishment of 

6 larch regeneration and reduces between species competition. Mineralized burned surfaces, enriched with 

7 nutrients are favorable for gemination of the small and light weight larch seeds. The expanse of larch 

8 forests is at present considered to be a carbon sink (Shvidenko et al 2007). However, an increase in fire 

9 frequency in response to observed climate changes in the area may result in conversion of this area to a 

10 source for greenhouse gases (1PCC 2007). Changes in air temperature, permafrost depth and extent may 

1 1 affect wildfire frequency (Kharuk et al, 2008). In spite of the fact that larch dominated forests occupy 

12 about 70% of permafrost areas in Siberia, data on fire occurrence in larch forests is presented by only a 

13 few publications (Vaganov and Arbatskaya 1996, Kovacs et al 2004; Kharuk et al 2005, 2008; 

14 Schepaschenko et al, 2008; Wallenius et al, 2011). For larch dominated areas of Central Siberia average 

15 FR1 was found 82±7 years (~64°N latitude), and 77±20 years for the southward “larch-mixed taiga” 

16 ecotone. For the northern boundary of larch forests (~71°N) FR1 value was estimated to be 320 ±50 years 

17 (Kharuk et al 2011). For the north-east larch forests of Siberia FRI was found to be, depending on site, 

18 50-80 and 80-120 yr (Schepaschenko et al, 2008). For southern Central Siberia Wallenius et al. (2011) 

19 reported a gradual increase of FRI from 52 years in the 18th century to 164 years in the 20th. 

20 The purpose of this work is to investigate wildfire occurrence in the central part of larch dominated 

21 communities (figure 1). 

22 2. Study area 

23 The test sites were located within the Kochechum River watershed. The Kochechum River is a tributary 

24 of the Nizhnyaya (lower) Tunguska River which turns northwest and flows into the Yenesei River. This is 

25 area is the northern part of the central Siberian plateau with gentle hills with elevations up to 1 000 m 

26 (figure 1 ). 
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This is a permafrost area with a severe continental climate. Mean summer (JJA) and winter (DJF) 
temperatures are +12°C and -35°C respectively. Mean summer and annual precipitation are 180 and 390 
mm, respectively (these values are means over the years 1997-2006 averaged for 1.5° x 2.5° grid cells and 
covering all test sites; figure 1; WWW1). 

The forests are composed of larch ( Larix gmelini Rupr.) with a mixture of birch ( Betula pendula 
Roth). Typical ground cover is composed of lichen and moss. Bushes were represented by Betula nana, 
Salix sp, Ribes sp, Rosa sp., Juniperus sp, Vaccinium sp, and Ledum palustre L. (Labrador tea). 

The wildfires across this landscape occur as ground fires due to low crown closure. The seasonal 
fires distribution is a single-mode (late May-June) with rare late i.e., (August-early September) fires 
(Soffonov et al 1999, Kharuk et al 2007). Periodic stand-replacing fires cause a mosaic of the (semi) 
even-age stands, embedded with older trees that survived the fire. 


60 ° 70 ° 80 ° 90 ° 100 ° 110 ° 120 ° 130 ° 



Figure 1 . Map of north central Siberia with the area of this investigation shown as a rectangle. Areas 
marked 1-3 are sites of earlier studies (Kharuk et al 2007, 2008, Kharuk et al 2011). Inset: Landsat image 


5 


1 showing locations (1-9) of test sites along the Kochechum River. Measurements of older trees from a 

2 study of long term wildfire trends were acquired in areas denoted by boxes. 

3 3. Materials and method 

4 3.1. Tree sampling 

5 The samples were collected in 200 1 and 2007 within the burned areas of larch forests along the river. 

6 There are no roads within the study area, and rivers provide the best access (figure 1). Temporary test 

7 sites were established within burned areas within about 1.0 km from the river and within a 160-420 m 

8 elevation range. Disks of larch bole cross sections for tree ring analysis were cut at the root neck level. 

9 Sampled trees included specimens of the dominant (even-age) trees, as well as older trees that survived 

10 stand replacing fires. 

1 1 Trees within test sites 1-9 (figure 1) were sampled and used for the fire return interval (FR1) 

12 analysis. In addition to this, older trees on supplementary sites were sampled (figure 1) for the purpose of 

13 estimation of long-term trends in fire frequency. The sample consisted of 58 trees. 

14 3.2. Sample analysis 

15 Tree ring widths were measured with a precision of 0.01 mm using the well-known LEMTAB-lll 

1 6 instrument. The dates of fires were estimated based on a master chronology constructed for northern 

17 larch forests (Naurzbaev et al 2004). Mean correlation with the master chronology was 0.54 which is 

18 satisfactory for our purposes. The COFECHA (Flolmes 1983) and TSAP (Rinn 1996) programs were used 

19 to detect double counted and missing rings. Fire-caused tree ring deletions were found only for 3 cases 

20 (out of about 75 analyzed). Relevant dates of fire events were adjusted based upon the master 

2 1 chronology. 

22 3.3. Fire return interval calculation 

23 Fire return intervals (FR1) were routinely estimated by tree ring calculation between consecutive fire 

24 scars: D, - D t _i, where D„ D U1 - dates of i and i-1 fires. Since many sampled trees have only a single bum 

25 mark, the dates of tree natality were included into the FR1 calculation where appropriate as described 

26 below. 
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1 It is known that within larch-dominated communities fires are mostly stand-replacing, and 

2 promote the formation of even-aged stands. Fresh bums with mineralized soil are quickly occupied by 

3 dense regeneration. Over time, development of a thick moss and lichen cover limits larch regeneration, 

4 (Kharuk et al 2008). Larch produce seeds annually with good harvests occurring on 5-7 year cycles 

5 (Forest Ecosystems 2002). Furthermore, 2-3 year old cones are known to be viable seed sources 

6 (Sofronov et al 1999). Importantly, ground fires regularly do not damage cones, leaving fire -killed stands 

7 as a source of seed. Thus both, tree natality dates and the dates of fires were used for FRI calculations. 

8 FRI were determined for each individual tree within test sites and then the mean FRI for a given site was 

9 calculated. 

10 Long-term history was based on fire events were for the period spanning the 17th century to 2007. To 

1 1 exclude the impact of a decreasing sample size for older stands i.e., “fading effect” the sample size was 

12 adjusted by tree age within groups. Stand replacing and non-stand replacing fires were investigated. 

13 

14 Along the burnmarks, “tree ring width growth accelerations” were also determined. “Growth 

1 5 accelerations” were identified by the following procedure. 1 . Expert visual analysis of increase in tree 

16 ring width.2. Calculation of the mean tree ring width 20 years before and 20 years after the increases 

1 7 began. This reference period (20 yr) approximately corresponds to the period of post-fire tree ring 

1 8 increment increases. 3. If the ratio (mean tree ring width after begining of tree width increase/mean tree 

19 ring width before tree width increase) > 2.0, the observed “growth acceleration” was considered 

20 significant. 4The date of the “growth increase” was compared with the air temperature record. It is known 

2 1 that for the boreal-forest zone, radial growth has been closely connected with temperature variability 

22 (Esper et al., 2010) . If that date coincided with air temperature increase, that particular “growth 

23 increment increase” was excluded from furthur analysis (because of the possibility of climate-driven tree 

24 ring increase). The above-mentioned “growth acceleration” dates were considered as a possible 

25 additional indirect sign of wildfire impact. It’s known that trees that survive fires (including those which 

26 do not have burnmarks) experienced a period of growth increase (e.g., Sofronov et al, 1999). 

27 
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1 Error analysis 

2 Including the tree natality date into FR1 analysis increased errors due to the 0-5 years lag of the post-fire 

3 tree establishment. The natality date also has to be adjusted to the “stump age”, i. e., the difference of real 

4 and measured stump height tree age. Even if a tree was cut at the root neck level, this difference can be 2— 

5 5 yr. In some cases disks sampled at the root neck level were not suitable for analysis; in these cases disks 

6 were sampled higher up the bole. This procedure entailed about 5 additional years of uncertainty. The 

7 post-fire regeneration had a high growth increment the first 1 5-20 yr, which gradually decreased due to 

8 increasing competition and decreasing active root layer depth (Soffonov et al 1999). In summary, the 

9 maximum total error was estimated to be 15 years. 

10 4. Results 

11 4.1. Fire return intervals 

12 Dates of tree natality and fire events are presented in figure 2. FR1 values for different test sites are 

13 considerably variable, i.e, from 131 to 349 years (with mean of 200±51 year; table 1). 
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Figure 2. Fire events chronology for the test sites 1-9 (locations are shown in figure 1). Light bars 
indicate tree natality dates, heavy bars indicate fire events, and diamonds indicate dates of initiation of 
radial growth acceleration. 
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Table 1 . Mean FR1 for nine study sites. 


Test sites # 

Mean FR1 

Sample size (N) 

1 

134 

5 

2 

203 

2 

3 

138 

9 

4 

165 

2 

5 

349 

5 

6 

278 

17 

7 

142 

4 

8 

256 

3 

9 

131 

4 

Overall Mean 

200±5 1 (p>0.05) 



4.2. Long-term fire events history 

The long-term fire event history (for the 17th through 20th centuries) was developed based on data for 
trees with natality dates before 1800, 1700 and 1600 yr, respectively (figure 3, table 2). To exclude the 
fading effect data were adjusted for tree natality dates: >200 yr, >300 yr and >400yr for 19th to 20th, 18th 
to 20th, and 1 7th to 20th century comparisons, respectively. 
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The number of wildfires in the 20th century increased relative to the 19th century (13 vs 8) (table 
2; N=58). For 18th, 19th and 20th centuries (N=33) fires number were 1, 6 and 9, respectively (table 2; 
N=33). Results for period 17th-20th centuries also showed an increase of fire events from the 17th to 
20th centuries (table 2; N=14). The minimum number of wildfires coincided with the air temperature 
decrease during the Little lee Age period (i.e., early 17th to early 19th centuries; figure 4). 
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Figure 4. Individual (gray lines) and running mean (with 50 yr window; dense solid line) radial increment 
data of sampled trees (N=58), dates of fire events (upper scale), and air temperature deviations for the 
northern Siberia and northern hemisphere (grey thin and bold solid lines; Briffa 2000). 


Table 2. Number of wildfires during 17 th - 20 th centuries 


Century 


Number of wildfires during 


Air temperature deviations, °C 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


12 


20th 

13 

a given century 
9 

6 

19th 

8 

6 

2 

18th 

- 

1 

0 

17th 

- 

- 

2 

Sample size (N) 

58 

33 
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5. Discussion 

5.7. FR1 

The majority of the fires in the larch communities we studied are stand-replacing. This was revealed by 
the fine-grained mosaic of (semi) even age stands within forested landscapes. For the types of forest 
communities the mean FR1 can only be reconstructed using age -class analysis or fire records (Sannikov 
and Goldammer 1996, Johnson and Miyanishi 2001). However, in our case the presence of fire scars 
indicated that non-stand replacing fires also occurred (figure 2). 

FR1 data for different sites within the study area (figure 1 ) differed by more than two times, from 
13 1 to 349 years with a mean of 200±5 1 years (table 1). The reason for this is that fires are rare events 
within the investigated area. The other cause is the non-uniformity of the topography and land cover 
within the study area. It is known that fire frequency is different for sunlit and shadowed slopes, as well 
as for bog areas (Beaty and Taylor 2001, Rollins et al 2002; Kharuk et al 2005, 2008). Meanwhile, even 
the lowest FR1 value (131 years) considerably exceeds reported data for adjacent southward forests (80- 
90 years; Vaganov and Arbatskaya 1996, Kharuk et al 2005). Mean FR1 values (200±51yr) exceeds 
published FR1 values for the boreal conifer forests (60-150 years: Payette 1992, Larsen 1997, Swetnam 
1996, Sannikov and Goldammer 1996). Very long FRls (up to 300 years) were reported for fire-protected 
forests in Europe and North America (Weir et al 2000, Heyerdahl et al 2001, Bergeron et al 2004, 
Buechling and Baker 2004). Evidently this is not the case, because within our study area fires were never 
suppressed. Low fire frequency is not favorable for the larch forests, because fires promote larch 
regeneration growth, i.e., larch is a “pyrophytic” species. The main tree growth constraints are 
permafrost thawing depth and soil drainage. Depth of seasonal thawing is dependent on exposition, moss- 
and-lichen layer thickness, and fire history. Fires not only increase permafrost thawing depth but also 
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increase soil drainage which is very important to larch growth. With time, an increase in the thermal 
insulator layer composed of the on-ground moss and lichen cover caused upward migration of the 
permafrost layer, and compression of the active root zone within a progressively decreasing upper layer. 
Fires also thin regeneration, decreasing within-species competition, and, thus, promoting tree growth 
because larch is an extremely shade-intolerant species. 

5.2. Long-term trends in wildfire number 

The advanced age of some sampled trees (>400 years; figure 3) allows estimation of the fire history 
within the study area back to the Little lee Age period. Dendrochronology data showed that cooling 
during the Little Ice Age period caused depression in the tree annual radial growth (figure 4). Fire event 
numbers since the 17th century approximately coincided with air temperature deviations, increasing with 
warming since the second half of 19th century (figure 4, table 2). For example, fires numbers increased 
from 8 in the 1 9th century to 1 3 in the 20th century. This phenomenon could not be attributed to 
decreased samples of older trees, i.e., “fading effect”, since sample sizes were adjusted by the natality 
dates. 

Local asynchrony of radial increment growth and temperature deviations on Fig. 4 could be 
attributed to the fire -induced increase of radial growth, as well as to the increment decrease during lag 
between fire event and increment increase begin (which is about 7-10 yr, Fig. 5). Thus, in the middle of 
20th century wildfires were observed within the majority of the test sites (Fig. 2). The other reason is a 
“divergence phenomenon”, i.e. growth-vs. -temperature divergence (D'Arrigo et al, 2008; Esper et al, 
2010 ). 

These estimates coincide with earlier reported data on fire frequency increase in the 20th vs 19th 
centuries for the southerly larch and mixed forests (sites 2 and 1 on the figure 1 , respectively; Kharuk et al 
2008). The causes of this trend can be both, natural and anthropogenic. Earlier it was shown that for site 1 
(figure 1) the FR1 decrease was caused by both, natural (warming) and anthropogenic causes. On site 2 
(figure 1) FR1 decrease was attributed to temperature increase mainly because the leading factor of fire 
ignition on the remote northern forests is lightning (>90% of cases in the northern Evenkia, whereas 
within southern forests >80% of fires are anthropogenic -caused; Kovach et al 2004). This is also true for 
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our study site since the remoteness of the area, as well as low population density in general; most of the 
fires are of natural origin. Similar observations (FR1 increase since Little Ice Age) were made for the 
northern boundary of larch forests (site 3 on figure 1; Kharuk et al 2011). All these data support the 
suggestion that observed climate change will lead to an increase in fire frequency (Gillett et a l 2004, 
Bergeron et al 2004, Girardin et al 2009). Comparison of FR1 along the meridian (figure 1) showed a 
northward increase of FR1: 77±20 years at ~61°N (site 1), 82±7 at 64°N (site 2), 200±51 years at about 
66°N, and 320±50 years at the northern boundary of larch dominant communities (71°N, site 3). The 
main reason of the northward FR1 increase is less incoming solar radiation and, consequently, shortening 
of the fire-danger period. 

5.3. Burns as a “simulator” of future warming 

Areas experiencing wildfires may be considered as a simulator of predicted warming impacts on northern 
forests. Indeed, in addition to soil enrichment with nutrients and decreased competition, fires cause an 
increase of permafrost thawing depth by a factor of 3 to 5 (Kharuk et al 2008), and increased soil 
drainage. Trees that survived wildfire showed an approximately double (1.93; p>0.95) increase of radial 
increment (figure 5). Comparison was made for the period for 25 (pO|[kjri]??) years before and after the 
fire (periods of +/-5 yr around zero were deleted to exclude effects of direct fire damage on trees. Some 
trees showed an extremely high response to fire affects. For example, the tree cross section shown on 
figure 5 showed an increase of radial growth about ten times in comparison with background 
observations. This tree was sampled at a latitude near the Polar Circle. 

Generally speaking, growth increases following fire scars should be measured along radii very 
most distant from the wound itself. But in our case we compare 37 specimens with the same 
pattern of fire damage; one of which shows outstanding increment growth (about ten times 
higher than background set; Fig. 5). The basic difference of this specimen with others sampled was the 
depth of soil thawing (about 1.5 m vs <0.3-0. 5 m for the other specimens), and good drainage since that 
tree was growing on the southern river bank. It is known that larch prefer drained soils (Schepaschenko 
et al, 2008). The observed radial growth increase (and, consequently, increased carbon sequestration) may 
be an alternative to the scenario of forests in climate-induced permafrost transformation areas becoming a 
greenhouse gases source (1PCC 2007). Thus, the vegetation dynamics and productivity of the burned 
areas deserves future investigations. 
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Figure 5. Diagrams for (1) individual tree ring widths (N=36) and (3) mean tree ring widths before and 
after fire; (2) - tree ring width for a specimen with an extremely high post-fire growth increment (inset). 
Data were compiled based on fires in 20th century (figure 2). Dates of fires were set to “zero” point. Note 
that post-fire growth increase has a lag about 7-10 yr. 

5.4. Indirect sign of wildfires 

Tree ring growth history showed periods of radial growth accelerations (i.e., tree ring width increase; 
figure 2-4). These increases are commonly considered to be climate driven (e.g., Shiyatov 2003) and may 
also contain information on fire events. The observed growth accelerations may also be caused by the 
above mentioned fire-caused soil melioration. Differentiation of these effects could be made based on the 
comparison with air temperature anomalies. The fire-induced origin of the acceleration is supported by 
the fact that in some cases the date of accelerations coincides with the bum marks on the trees from the 


same test site (figure 2).These effects deserved future investigations based on a larger sample sizes. 
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6. Conclusions 

FR1 on the study area is considerably longer than in southerly territories. Along the 100 degree meridian 
(figure 1) FR1 values increased northward, 77±20 years at~61°N (site 1), 82±7 years at 64°N (site 2), 
200±51 years at about 66°N, and 320 ±50 years at the northern boundary of larch stands (~71°N). The 
number of fire events during the Little lee Age period (17-1 8th centuries) was approximately half the 
number observed in 19-20th centuries. Fire -caused soil melioration (basically soil drainage and thawing 
depth increases) caused a radial growth increase about 2 times (with >6 times in extremes). This effect 
may simulate predicted warming impact on the larch growth in the permafrost zone. 
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